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ABSTRACT 

We present the first direct measurement of tlie spatial cross-correlation function of high-mass X-ray 
binaries (HMXBs) and active OB star-forming complexes in the Milky Way. This result relied on a 
sample containing 79 hard X-ray selected HMXBs and 458 OB associations. Clustering between the two 
populations is detected with a significance above 7cr for distances <1 kpc. Thus, HMXBs closely trace 
the underlying distribution of the massive star-forming regions that are expected to produce the progenitor 
stars of HMXBs. The average offset of 0.4+0.2 kpc between HMXBs and OB associations is consistent 
with being due to natal kicks at velocities of the order of 100+50kms"'. The characteristic scale of the 
correlation function suggests an average kinematical age (since the supernova phase) of ~4 Myr for the 
HMXB population. Despite being derived from a global view of our Galaxy, these signatures of HMXB 
evolution are consistent with theoretical expectations as well as observations of individual objects. 

Subject headings: Galaxy: evolution, open clusters and associations, stellar content, structure ; Stars; emission-line. 
Be, neutron, supergiants ; X-rays: binaries 



1. Introduction 

High-mass X-ray binaries (HMXBs) are systems 
in which a compact object (usually a neutron star, but 
sometimes a black hole candidate) accretes from a 
massive stellar companion (M > 10 Mq). The mass- 
age relation in stellar evolution predicts that ~ 10^ yr 
elaps es between starbirth a nd supernova in high-mass 
stars dSchaller et al.lll992l) . Thus, HMXBs are rela- 



tively young systems which are not expected to mi- 
grate far from their birthplaces: sites with a recent 



history of massive star formation, i.e., the OB as- 
sociations that trace the Galactic spiral arms (e.g. 
.Brown et al. 1999,) . 

Observational evidence linking Galactic HMXBs 
and OB associations has been demonstrated in indi- 
vidual cases in which the connection is attested by 
other factors: i.e., consistent proper motions and dis- 
tances, the position of the HMXB donor star with re- 
spect to the main sequence (MS) of the OB associa- 
tion on color-magnitude diagrams, etc. On a Galactic 
scale, however, the evidence is limited to comparing 
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distributions of Galact ic longitudes (/) or galactocen 
trie d i stanc es (R) (e.g., 'Grimm et al. ^1 20021: iDean et al, 



20051: iLutovinQ y et al 2005; Boda ghee et al.l l2007t 



Lutovinov et al. 2008). Both HMXBs and OB as- 



sociations show maxima around Galactic longitude 
; ~ 30° which corresponds roughly to the direction 
of the Norma and Inner Perseus Arms (upper panel 
of Fig.[TJ. Moderate peaks are found towards the Ca- 
rina, Scutum, and Sagittarius Arms. This suggests 
that HMXBs are spatially correlated with active sites 
of massive-star formation. This correlation is further 
exemplified by the similarity in the distributions of 
galactocentric distances of HMXBs and OB associa- 
tions: both distributions can be approximated by Gaus- 
sians centered on the galactocentric distance of the Sun 
(lower panel of Fig. [Til. 

Recent efforts to describe the Galactic distribu- 
tion of HMXBs raised the intriguing possibility of 
an offset between the peak in the HMXB longitu- 
dinal distrib ution with respect to the tangent to the 
Norma Arm (Lutovinov et al. 20051: Dean et al.ll2005 



Bodaghee et al.lbOOTtlLutovinov et al.ll2008l) . This di 
rection features the highest formation rate of massive 
stars (Bronfman et al. 2000, ). An offset would im- 
ply that we are witnessing the systematic delay (= 
kinematical age) between the epoch of star formation 
and the moment when these massive binaries became 
X-ray emitters. Figure|2] presents the population of 
HMXBs and OB associations in the Milky Way as it 
would appear to an observer situated above the Galac- 
tic Plane. Thanks to an increase in HMXB statistics. 



and by using a recent Galactic spiral arm model dVallee 
I2OO8), we can now attribute the apparent offset, at least 
partly, to foreground sources that are unrelated to the 
Norma Arm, and to sources located towards the tan- 
gent to the Inner Perseus Arm. The extension of the 
latter arm into the Inner Galaxy was not considered i n 
previous spiral arm models (see also Hou et all2009l) . 

Clearly, a precise study of the spatial relation of 
massive star-forming regions and the HMXBs they 
spawn can offer valuable insight into stellar and Galac- 
tic evolution. The drawback with longitudinal or 
galactocentric distance distributions is that these his- 
tograms consider only a single dimension at a time. In 
effect, they project an entire spatial component {R in 
the case of a longitudinal distribution, and / in the case 
of a galactocentric distance distribution) onto the other 
component. Precious information about the true prox- 
imity of objects in the two populations is lost. Hence, 
these methods can not accurately assess the spatial 
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Fig. 1. — Distribution of Galactic longitudes and 
galactocentric distances (when known) of HMXBs 
(shade d histogram) and OB associations from lRusseil 
( l2003b (thick curve: divided by 5). The locations of 
the spiral arm tangents are shown. The distributions of 
HMXBs in these projected spatial directions are com- 
patible with those of OB associations. 



relation, restricting their ability to provide the char- 
acteristic scales that would help us understand how 
HMXBs are distributed in the Milky Way. 

A better approach is to consider the spatial (or 
cross-) correlation function ^ (or SCF). This function 
provides a statistical measure of the clustering of point 
sources in multi-dimensional maps of the sky. It is fre- 
quently used to relate active galactic nuclei with clus- 
ters o f galaxies (e.g. [Peebles 1980; Landy & Szalavl 
1993b . In Section[3] we adapt the methods described in 
these papers to suit the study of HMXBs on a Galactic 
scale, and we assess biases in SectionH] Our findings 
are discussed in the context of stellar and Galactic evo- 
lution in Section|5j 



2. Data 



We selected all sources from lBodaghee et alJ(l2007h 
which are confirmed or strongly suspected of being 
HMXBs. This catalog contains hard X-ray detected 
sources (>20keV), meaning that their detection does 
not depend on the photoelectric absorption that blocks 
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Fig. 2. — Galacti c distribution o f HMXBs whose distances are known (79, filled triangles) and the locations of OB 
associations from lRusseill ( l2003b (458, circles). The symbol size of the latter is proportional to the amount of activity 
in the associatio n (= amount of ionizing photons as determined from the radio continuum flux). The spiral arm model 
of IValled (l2008h is overlaid with the Sun situated at 7.6 kpc from the Galactic Center (GC). HMXBs whose d i stance s 
are not known have been placed at 7.6 kpc (23, empty triangles), i.e., the Sun-GC distance assumed in I Vailed (120081) : 
note that these sources are not included in the cross-correlation analysis. The shaded histogram represents the number 
of HMXBs in each 15° bin of galactic longitude as viewed from the Sun. 



soft X-rays (<5keV). A number of discoveries have 
been made since the original catalog was published, 
and these new sources were included as well. Many 
HMXBs are listed on the IGR Sources Page which we 



keep currenfl The date of April 1, 2011, represents 
the cut-off after which we no longer considered new 



'http : //irfu . cea . f r/Sap/IGR-Sources 
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information. 

For each HMXB, we searched through the literature 
and extracted the most recent distance measurement 
and related uncertainty when available. We strove to 
consider only distances that were derived from ob- 
servations of the optical/IR counterpart, and not dis- 
tances assumed from, e.g., being in the line of sight 
of a particular spiral arm or the Galactic Center (GC). 
Any HMXB suspected of or confirmed to belong to 
the Magellanic Clouds was excluded. The final sam- 
ple consists of the 79 X-ray selected HMXBs listed 
in Table[T] The positions a nd distances of OB-forming 
complexes are drawn from Russeill ( 2003 ). These com- 
plexes are large and excited H II regions which include 
OB associations, but we will use the term "complex" 
interchangeably with "association." We rejected OB 
associations that are located more than 20kpc from the 
GC resulting in 458 associations. 

We emphasize that the HMXB and O B samples are 



drawn exclusively from the Milky Way. iLinden et al. 

argue that galaxies with a metallicity compa- 
rable to that of the Milky Way produce few luminous 
blac k hole HMXBs suc h as CygX-1 or SS433 (see 
also iMirabel et all 1201 lb . Hence, the HMXB sample 
contains mostly neutron star HMXBs classified into 
one of three groups: the transient emitters with Be 
companions (BEXBs), the mostly persistent emitters 
with supergiant companions (SGXBs), and the fast X- 
ray transients with supergiant companions (SFXTs). 

3. The Spatial Correlation Function 

For a given HMXB in a volume element 6V], the 
probability 6P (in excess of Poisson) of finding an OB 
association in a volume element SV2 separated by a 
distance r (in kpc) is: 



6P = nin2[l +^{r)]6Vi5V2 



(1) 



The mean number densities are given by «i and 
n2 where the subscripts refer to the two populations: 
HMXBs (1) and OB associations (2). These are young 
objects that remain within 45 + 20 pc of the Galactic 
Plane jReed l2000l) . so we assume that the displace- 
ment along the z-axis is negligible, and surface (rather 
than volume) elements are used. 

Rings of 1 -kpc thickness (from to 20 kpc) are con- 
structed around each HMXB. In each ring, we count 
the number of observed OB associations, and we keep 
a separate tally of the number of OB associations 
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Fig. 3. — Spatial correlation functions ^(r), where r 
represents the distance from a given HMXB (in kpc), 
are presented for both estimators. Three random dis- 
tributions were considered: a Gaussian ring centered 
at 7.6 kpc with cr. - 2.5 kpc; a uniform disk for 
< r < 16 kpc; and a uniform disk with a hole at 2 kpc 
from the GC. The curves are slightly shifted horizon- 
tally in order to allow the 1-cr error bars to be clearly 
distinguished. At short distances from a given HMXB, 
^ > 1. This suggests that the closest neighbors of 
an HMXB tend to be observed OB associations rather 
than OB associations drawn from a random distribu- 
tion. 

drawn from a random distribution. An HMXB-OB 
pair is referred to as a DD^-pair (for data-data pair), 
while an HMXB-random OB pair is labelled DR12 (for 
data-random). In this way, we co nstruct f for eac h ra- 
dius according to the definition of lPeebles (Il980h : 



HRDDn 
syn 1 

noDRn 



(2) 



If ^ = 0, which implies that each ring contains as 
many DZJ-pairs as D/? -pairs, then Eq.[T]is simply a uni- 
form Poissonian probability. However, if ^ > 0, then 
there is a higher chance of an HMXB having a neigh- 
bor that is an observed (rather than a randomized) OB 



association. 



Landv & Szalavl (Il993h propose a more 



robust estimator for ^ in which the variance is nearly 
Poissonian: 
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Fig. 4. — Distributions of distances between pairs 
of HMXBs and observed OB associations (red 
shaded histogram), and between pairs of HMXBs 
and Gaussian-randomized OB associations (divided by 
10^, blue histogram). A KS-test (lower panel) yields a 
probability < 10"^ of statistical compatibility. 



DDx2 - DRn - DRi\ + RRxi 

RR]'y 



(3) 



Three random distributions were considered: a disk 
model in which the OB associations were uniformly 
distributed between an 16 kpc from the GC; a disk 
model with a hole, i.e., a uniform distribution from 
2 kpc to 16 kpc from the GC; and a Gaussian distri- 
bution centered at 7.6 kpc from the GC with cr. - 
2.5 kpc. The boundaries, centroid, and width of the 
Gaussian were chosen to be consistent with the ob- 
served distribution of HMXBs and OB associations 
(see Fig. [Til. Each random distribution contained the 
same number of OB associations as the observed OB 
sample (458) occupying similar surface areas so no 
and n« are equal. One thousand trials were performed 
with each randomization model. The average ^, along 
with its 1-cr uncertainties, are shown in Fig. [3] 

All three randomization models and both estimators 
of ^ lead to the same result: near an HMXB, the prob- 
ability of finding a known OB association is higher 
than expected from Poisson statistics. In other words, 
HMXBs and OB associations are clustered together 



Our measured value of ^(r < 1 kpc) = 1.61 + 0.22 
(Gaussian case) implies that xDD ~ DR (where x ~ 
2.6), i.e., around a given HMXB, one counts between 2 
and 3 times as many observed OB associations as ran- 
domized ones. We are well below the value of ~458 
anticipated for perfectly-correlated samples, but statis- 
tically well above the value expected for no correlation 
whatsoever (^ ~ 0). 



For r < 3 kpc, the clustering signal from the lPeebles 
esti mator is 7-1 Icr in excess of Poisson (up to 
17cr for the iLandv & Szalavl (119931) estimator). The 



observed and random surface-density distributions are 
statistically compatible only at large radii (i.e. r > 
5 kpc). Notice that while ^{r < 1 kpc) > ^{2 < r < 
3 kpc), the standard deviation of ^ decreases as we 
move away from an HMXB (up to a certain limit). 
This is simply a consequence of the increase in the 
surface areas, and hence the number of possible pairs, 
with increasing radius. 

Of the three models, the Gaussian distribution pro- 
vides the best approximation since its ^ deviates the 
least from Poisson statistics at short distances. The 
disk distributions (with and without a central hole) are 
no longer considered. 

4. Statistical and Systematic Biases 

Observational biases affect the correlation function 
of the objects in the survey. This is not trivial as the 
distribution of HMXBs and OB associations are, a pri- 
ori, not well understood. The type of random dis- 
tribution that one chooses leads to a systematic bias. 
This is readily apparent from Fig.|4] where the dis- 
tances between pairs of HMXBs and observed OB as- 
sociations are compared with the distances between 
pairs of HMXBs and Gaussian-randomized OB associ- 
ations. Despite the fact that the distributions have sim- 
ilar mean and variance values, the higher moments of 
the distributions (skewness and kurtosis) disagree. A 
Kolmogorov-Smirnov test yields a probability of less 
than 10"^ of statistical compatibility between the two 
distributions. 

Selection bias is also important. An HMXB sit- 
uated far from the Sun needs to be correspondingly 
more luminous in order to be detected in the X-rays. 
Large distances also make the optical classification as 
an HMXB more difficult: line spectroscopy of the 
donor star is hindered by reddening and absorption 
from interstellar material in the Galactic Disk. This 
leads to a preponderance of HMXBs (as well as ob- 
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Fig. 5. — Same as Fig. |3]but restricted to objects within 
8 kpc from the Sun. The random catalog assumes 
a Gaussian-ring distribution. Clustering between 54 
HMXBs and 361 OB associations persists with a sta- 
tistical significance of 7cr. 

served OB associations) close to the Sun, and a paucity 
of such objects situated behind the GC and its bar: a re- 
gion referred to as the "zona Galactica incognita" (e.g.. 
Vailed I2OO2I) . The random distributions of HMXBs 
and OB associations that we generated did not con- 
sider this observational sampUng bias. 

It is reasonable to expect symmetry in the distribu- 
tion of HMXBs and OB associations around the Galac- 
tic Disk. However, the fact that there are less HMXBs 
detected (with distances measured) in the zona incog- 
nita will surely affect the SCF Therefore, we gener- 
ated the SCF for a restricted boundary corresponding 
to a circle of 8-kpc radius around the Sun — a perime- 
ter inside of which most HMXBs and OB associations 
should be detectable and their distances known with 
reasonable accuracy. When only pair counts of objects 
(54 HMXBs and 361 OB associations) in the Solar 
neighborhood are considered, the clustering signal at 
small radii persists with a statistical significance of 7cr 
for r < 1 kpc (see Fig.|5]l. 

Some HMXBs have distances with uncertainties 
that are on the order of, or larger than, the 1 -kpc spa- 
tial binning that we used. The counts per ring will 



Fig. 6. — The clustering remains significant even when 
the HMXBs are perturbed within a range defined by 
the uncertainties on their distances. Top: Three-sigma 
boundaries (gray shaded region) from multiple SCFs 
generated by randomly shuffling (according to a Gaus- 
sian distribution) the line-of-sight distance to each 
HMXB within its distance uncertainty. The blue dotted 
line and data points represent the non-perturbed (i.e., 
observed) SCF as presented in Fig. [3] Bottom: The dis- 
tribution of the clustering significance from all reshuf- 
fling trials where significance is defined as ^{r)/o-. 
From left to right, the panels present the distributions 
for f(r < 1 kpc), ^(1 < r < 2), and ^(2 < r < 3). 



fluctuate depending on the distance precision quoted 
in the literature. To account for this bias, we shuffled 
each HMXB according to a Gaussian profile within its 
line-of-sight distance range as defined by its error bars. 
When an error bar was not provided in the literature, 
we set the uncertainty range to be equal to the aver- 
age of the known error bars of the sample (1.5 kpc). 
Then, we regenerated ^ for each of the 100 randomly- 
shuffled HMXB distributions. Keep in mind that the 
OB distribution is itself randomized 10^ times per trial, 
which means that the total number of different HMXB- 
OB configurations that we tested is actually 10^. 

Even when the HMXBs are perturbed within their 
line-of-sight uncertainties, the clustering signal re- 
mains above 7 cr in all trials for r < 3 kpc from a 
given HMXB (Fig.|6]l. Consider a conservative as- 
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Fig. 7. — The SCF is presented comparing the lo- 
cations of OB associations, and an equivalent num- 
be r of points rep resenting the four-arm spiral model 
of lRusseill ( l2003b . The randomized catalog assumes a 
Gaussian ring distribution and the Galactocentric dis- 
tanceof the Su n is set to 8.5 kpc, i.e., the value used in 
Russeill (l2003b . 



sumption that says the clustering is significant only 
when ^ > 1 (i.e., when each ring around an HMXB 
contains at least twice as many observed OB associa- 
tions as random OB associations). We find that 97% 
of the correlation functions describing the shuffled dis- 
tributions satisfy this stringent criterion. A little over 
half of the trials (53%) yield a ^ value that is above this 
cutoff at the 3-cr level. Based on the >99% confidence 
interval displayed in Fig.|6] we estimate a systematic 
uncertainty of 5% on ^. 

Potential differences in the scale heights of the two 
populations has thus far been ignored. A few HMXBs 
present high velocitie s tangential to the Galacti c Plane 



(e.g., 4U 1907+09, iGvaramadze et alJ 120111) . By 



projecting this velocity component onto the Galactic 
Plane, we can expect to underestimate the true offset 
between an HMXB and a nearby OB association. This 
will lead to samples that are more correlated than they 
should be. Therefore, we calculated the Galactic scale 
height for members of our HMXB and OB samples. 
The average scale heights are in good agreement at 
0. 1 kpc for both samples with standard deviations of 



Fig. 8. — The SCF is presented for the case in which 
a randomized (Gaussian ring) HMXB catalog is taken 
as the observed set and is then compared with the OB 
associations. 



0.1 and 0.4 for HMXBs and OBs, respectively. This 
suggests that at the 3-cr level, HMXBs are located at 
most 0.4 kpc above or below the Plane. This value 
is consistent with the characteristic scale and average 
minimum offset that are discussed in Section|5] Ninety 
percent of the HMXBs (70/79) are within 0.25 kpc 
of the Plane with no HMXBs located further than 
0.55 kpc from the Plane. Therefore, we are justified 
in using the 2-D distribution since any increase in the 
correlation significance caused by this projection ef- 
fect will be negligible given that these offsets are less 
than the 1-kpc width of each ring in our correlation 
analysis. 

We tested the correlation between HMXBs and OB 



associations against the spiral arm model of iRusseil 
(l2003b whose equations represent the best-fitting four- 
arm logarithmic spiral to the locations of the star- 
forming complexes shown as blue circles in Fig.|2] 
Here, we adopted the Sun-GC distance of 8.5 kpc as- 
sumed in the model. This model yields an overdensity 
of points describing an arm in the inner Galaxy com- 
pared with the same arm in the outer Galaxy (i.e., as 
one loops through the range of angles that trace an arm 
from the inner to the outer Galaxy, the actual distance 
separating successive points becomes larger). Thus, in 
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Fig. 9.— The SCF compares HMXB s with 133 glob- 
ular clusters from iBica et alj (120061) that are located 
within 20 kpc of the GC. The random catalog of glob- 
ular clusters assumes an exponential decay law. The 
SCF is consistent with Poisson statistics at close radii, 
which means that the neighbors of HMXBs are just as 
likely to be drawn from the observed distribution as 
they are from a random one: i.e., there is no clustering 
between HMXBs and globular clusters. 

order to reduce this bias, we extracted 115 points cho- 
sen at random along each of the four arms (460 total: 
i.e. equivalent in number to the 458 observed OB as- 
sociations). We performed a visual check and verified 
that the spiral pattern was easily discernible with no 
major gaps (> 2 kpc) along an arm. 

Figure]?] shows that there is no significant correla- 
tion between OB associations and the spiral arms with 
^(r) ~ for all r (i.e., DD ~ DR for all r). Similarly, 
we found no significant clustering between HMXBs 
and the points representing the spiral arms. Restrict- 
ing the analysis to objects in the Solar vicinity (i.e., 
within 8 kpc of the Sun) leads to similar conclusions, 
so we can rule out an apparent lack of distant OB as- 
sociations (or HMXBs) skewing the results. In other 
words, around any given OB association (or HMXB), 
one tends to find as many neighbors drawn from a 
randomized Gaussian ring distribution, as neighbors 
representing the spiral arms. These comparisons of 
our observed stellar samples with the Galactic mod- 



els all ow us to conclude, as acknowledged by iRusseil 
(120031) . that such spiral arm models are an overly sim- 
plistic representation of the real picture. 

In another test, we randomized the HMXB distri- 
bution (using the Gaussian ring approximation) and 
we assumed that this sample was the observed HMXB 
set. The ^ relating randomized HMXBs and OB as- 
sociations is presented in Fig. [8] There is a moder- 
ate deviation from Poisson statistics at short radii for 
the Peebles (Il980l) estimator (whose values are >0 at 
all radii, as expected since the distributions are simi- 
lar), w hile there is no deviation for the Landy & Szalayl 
(Il993h estimator. This suggests that even though the 
randomization model mimics the global distribution of 
the observed HMXB sample, it is the observed HMXB 
distribution that is clustered with the OB associations. 

Clustering between HMXBs and OB associations is 
expected. How would ^ react if HMXBs were com- 
pared to a population of sources for which no spa- 
tial correlation is expected? We tested the correlation 
functions of HMXBs agains t a set of 133 globular clus- 
ters from Bica et al. I (l2006h that are located less than 
20 kpc from the GC. Globular clusters contain older 
populations such as cool KM dwarf stars and low- 
mass X-ray binaries (LMXBs). Unlike HMXBs, glob- 
ular clusters (and LMXBs) are densely packed in the 
Galactic Bulge and their numbers dr op exponentially 



with increasing radius from the GC dBodaghee et al. 



2007n . The random distribution of globular clusters 



was thus modeled as an exponential decay law ad- 
justed to be consistent with the observed distribution. 
Figure|9] shows that ^ between HMXBs and globular 
clusters is consistent with for r < 1 kpc from an 
HMXB. For any given HMXB, its immediate neigh- 
bors were just as likely to be drawn from the observed 
set as from the random set. Shot noise dominates at 
large distances. 

The results of these tests suggest that the cross- 
correlation function is assessing the real spatial clus- 
tering of populations (within the uncertainties) rather 
than some global property that they share. In other 
words, the significant clustering that we measured be- 
tween HMXBs and OB associations (Fig.O is due es- 
sentially to the physical locations of the individual ob- 
jects in space, as opposed to being an attribute of the 
overall populations, e.g., that they appear to trace the 
spiral arms, or that they seem to follow a Gaussian dis- 
tribution, etc. 



8 



Migration Distance [kpc] 




1111,1111 
/ 


, , , , ,1111, 

- 

\ 

\ 




\ 




\ 


J. - / 


















S \ 




\ \ 


^Cf<l kpc) 




__f(15:r<3) 




f(2Sr<3) 




V = SOOkms-i 




-5 5 

Migration Timescale [Myr] 



Fig. 10. — The values of the SCF that result from shift- 
ing the positions of the HMXBs along their Galactic 
orbits. The shifts correspond to 1-10 Myr (in steps of 
IMyr) of circular motion at a velocity of 200 km s ' 
in the forward (+) and reverse (-) directions of time. 
The continuous curve represents the values for ^(r < 
1 kpc), the dashed curve is ^(1 <r< 2), and the dotted 
curve is ^(2 < r < 3). 

5. Discussion 

There are several questions that we seek to address 
using the SCF. What is the characteristic scale of the 
clustering between HMXBs and OB associations? Can 
this scale be used to constrain the amount of migra- 
tion due to the effects of Galactic rotation (see, e.g., 
Ramachandran & Deshpanddl 19941) . and perhaps due 
to the primordial "kick" velocity imparted to the sys- 
tem during the formation of the compact object? Over- 
all, what can this scale tell us about the evolutionary 
history of HMXBs? 



5.1. Galactic Rotation 

Accor ding to the derisity-wave theor y of spiral 



galaxies dUn & Shulll964l:lLin et alJI 19691) . the veloc 



ity of the spiral pattern is only half that of the material 
moving in and out of the arms. These arms appear sta- 
tionary to an observer in a non-inertial reference frame 
that is rotating at the pattern speed, while the motions 
of individual objects within the Galaxy are not. As 



material encounters the quasistatic density waves, it is 
shocked and compressed which allows the Jeans crite- 
rion to be satisfied locally, triggering fragmentation of 
the cloud, and ultimately star formation. 

Massive stars in binary systems hav e a main se- 
quenc e (MS) lifetime of around 5x 10^ yr dSchaller et al 



19921) . The exact age depends strongly on the stellar 
mass and spectral type. In addition, there is a system- 
atic delay of ~ 4 x lO*" yr betw een the supernova phas e 
and the X-ray emission phase ( van den Heuvel 1983h . 
The X-ray phase lasts only 10"* yr ( Iben et all 1995 ) 
which is negligible on these timescales. In 10 Myr, 
a combination of global (Galactic rotation) and local 
effects (e.g., supernova kicks, cluster outflows, and 
dynamical ejections) will cause newly-formed O and 
B stars to eventually migrate away from the dense 
gaseous regions in which they were born. In principle. 



each spiral arm includes three components (iLin et al. 
1969h : a dusty strip in which the gas concentration 



and OB star formation rate are highest; a lane of lumi- 
nous and recently-formed OB stars that have ionized 
large regions of gas (HIl); and a wider lane of post- 
MS, red, or dying supergiant stars with smaller Hll 
regions. The HMXBs should be found among the 
recently-formed OB stars and the post-MS stars. 

Outside of the inner Galaxy (i.e., galactocen- 
tric radius R > 2 kpc), objects orbit around the 
GC a t a velocity of ~200kms~' (e.g.. iMerrifield 



19921: iBr and & Blitzl 1 19931: iGlushkova et alJ Il999t 



Sofue et al., ,2009) . Thus, a star at the Solar galacto- 



centric radius (7.6kpc) will have moved 6.3 x 10'^ km 
(200 pc) in 1 Myr (defined as the migration timescale 
or t), assuming a circular orbit. In this manner, we 
rotated the HMXBs in space along their galactocen- 
tric orbits corresponding to t = 1-10 Myr (in steps of 
1 Myr) of circular motion in both the forward (+) and 
reverse (-) directions of time (clockwise and counter- 
clockwise, respectively, in Fig.|2l). FigurefTOl presents 
the values of the SCF resulting from these shifts for 
small radii (r < 3 kpc from an HMXB). When the 
radius is larger than ~5 kpc, ^ is compatible with Pois- 
sonian statistics, and so these curves are omitted for 
clarity. 

Certain shifts increase the amplitude of the clus- 
tering signal. Increases in the amplitude of ^ are di- 
rectly related to an increase in the number of observed 
OB associations contained within the specified radius 
around the shifted HMXBs. This is because the ran- 
domized OB distribution is axisymmetric with respect 
to the GC, and so the number of random OB associa- 
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tions inside the radius remains relatively constant un- 
der shifts in either direction. Given that these shifts 
are small (< 2kpc), their effect on ^ diminishes with 
increasing distance from an HMXB, as illustrated by 
the trend towards flatter curves for ^(1 < r < 2kpc) 
and ^(2 < r < 3 kpc). We are mainly interested in 
clustering at short scales so, henceforth, our discus- 
sion will focus on ^(r < 1 kpc) (the solid curve in 
Fig.[T0ll. At smafl radii, the IPeeblesI (ll98Cll) estima- 
tor provides sufficiently reliable results while also be- 
ing more efficient comp utationally compared with the 
Landy & Szalayl (Il993b estimator. 

If HMXB migration were related solely to Galac- 
tic rotation, then we would expect to see the ampli- 
tude of ^ maximized for migration timescales around 
-lOMyr. Instead, we find that the amplitude of ^ 
is maximized for migration timescales between -2 
and -7 Myr. This suggests that young binary systems 
(prior to the formation of the compact object) are grav- 
itationally bound to their birthplaces, in most cases. 

Another consequence of Galactic rotation is that 
more HMXBs should appear on the leading edge of a 
spiral arm than on its trailing edge (because of the fac- 
tor 2 difference between the bulk and pattern speeds). 
Visual inspection of the Carina Arm (Fig.|2]l reveals 
no clear asymmetry in the distribution of HMXBs with 
equivalent numbers on either side of the arm. In the in- 
ner Galaxy, the tight spacing between arms combined 
with the uncertainty in the radial direction prevents us 
from visually assigning an HMXB to a specific arm. 
Nevertheless, since ^ is a good indicator of the average 
proximity of HMXBs to the OB associations that trace 
these arms, we should see an asymmetry in its distri- 
bution with higher values of ^ at negative timescales 
than at positive timescales. 

Indeed, there is a slight preference for migration to- 
wards the leading edge of spiral arms as hinted at by 
the asymmetric distribution and the downward trend 
in ^ for -5 Myr < r < +5 Myr. The amplitude of ^ is 
greater for negative migration timescales (r ~ -4 Myr) 
than it is at rest (r = OMyr). However, the ampli- 
tude of f for small positive migration timescales (i.e., 
+2 Myr < r < +3 Myr) is equivalent to that at rest 
which suggests that some HMXBs migrate in direc- 
tions different from that of the bulk motion. This is 
another indication that we are not witnessing the ef- 
fects of Galactic rotation alone. 



5.2. Kick Velocities 

Galactic rotation can not fully explain the be- 
havior of ^ under different migration timescales. 
One mechanism suspected of giving an HMXB a 
significant velocity is the primordial kick inherited 
by the binary during the supernova event that cre- 
ated the compact object. This runaway velocity 
(with respect to its OB associat ion) can be pro - 
duced by asymmetric supernovae (IShklovskii Il970l) . 
or from recoil due to anisotr opic mass los s from 
the primary to the secondary (lBlaauwlll96ll) . Fac- 
tors external to the binary can also allow it to es- 
cape from its parent associatio n: e.g,, dynamical 



ejection and cluster outflows ([Poveda et al.l 11967 
Pflamm-Altenburg & Kroupal201ol) 



Theoretical estimates fo r post-SN kick veloci- 
ties a re va ~70-120kms ' (iBrandt & Podsiadlowski 
1995b . assuming a non-eccentric, 1-25-d orbital pe- 



riod binary system composed of a 5 Mq star that ex- 
plodes in an asymmetric supernova leaving behind a 
1.4 Mq NS that tugs its hefty 15 Mq donor star Ve- 
locities typical of isolated pulsars (200 km s"' or even 
350 km s ' ) have been suggest ed depending on the ini- 
tial conditions in the binary ( Portegies ZwartI 19951: 



van Bever & Vanbeveren|[l997 ). van den Heuvel et al 



(120001) distinguish between kick velocities in SGXB 
systems from those of BEXBs: ~40 and ~15kms 
respectively. Selecting only the 38 SGXBs or only 
the 35 BEXBs from the HMXB sample yield consis- 
tent results: both populations show >6-cr clustering 
with OB associations with no statistically significant 
differences between SGXBs and BEXBs. 

In certain cases, the value of the kick velocity 
can be derived from the the radial velocity or proper 
motion when the identity of its parent OB associa- 
tion and the line-of-sight distances are known. The 
HMXB can then be traced back to its birthplace lead- 
ing to an estimation of its kinematical age. Obser- 
vations of the SGXB LS 5039 set v ~ 150 km s ' 
which constrain s the time sinc e its supernova to 
within 1.1 Myr (iRibo et al.1 [2002h . Another BEXB, 
LSI+61°303, has v = 27 km s~ ' and experienced it s 
supernova 1.7 + 0.7 Myr ago (iMirabel et al.l |2004 . 
The SGXB 4U 1700-377 shows v ^ 75kms-' wit h 
its supernova 2.0 + 0.5 Myr ago (lAnkav et al.ll200lb . 
VelaX-1, which is also a SGXB, has v = 90kms ' 
with 2 + 1 Myr having ela psed since its supernova 



(Ivan Rensbergen et al .119961) . [ Recently, the velocity of 
the SGXB 4U 1907 h-09 was measured at ~ 160 km s ' 
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with respect to t he Galactic Plane. Its kin ematical age 
is only 0.1 Myr (iGvaramadze et al.ll201 11) . The latter 
two systems present bow shocks as they move through 
the interstellar medium confirming the supersonic ve- 
locities of these HMXBs. Table|2] summarizes these 
measurements. 

We emphasize that there is no clearly-identified OB 
association for most HMXBs in our sample. Yet the 
vestiges of this dynamic history are imprinted on the 
overall HMXB distribution. Consider the distance sep- 
arati ng each HMXB fro m its nearest OB association 



(e.g.. Ivan den Heuvel et a l. 2000). Excluding separa- 



tions larger than 1 kpc, because they are less likely to 
be from post-natal kicks, we obtain an average min- 
imum distance of =0.4±0.2kpc. With this rmin, 
and with ^ maximized at t ~ -4 Myr (Fig.fTOli. we 
estimate an average velocity of 100±50kms"' for the 
HMXBs in our sample. This is consistent with mea- 
surements of the kick velocities in individual objects 



(see T able[2land, e.g., 
19791) . 



Gunn & Ostrike^ll970t [Stone 



Alternatively, with r^i„ - 0.4 + 0.2 kpc and as- 
suming V = 100 km s"', this translates to a migration 
timescale of 4+2 Myr which is consistent with theo- 
retical predictions for the average kinematical age of 
runaway massiv e binaries ( Blaauw 1961 ; Stone 19791; 
van Oiienll989l) . Notice that these kinematical age and 
distance scales are consistent with the distribution of 
^ whose value is maximized between -2 and -7 Myr 
(Fig.fTOb. The lack of an equivalent increase in the 
value of ^ between +2 and +7 Myr (which we would 
expect for kicks in random directions) is consistent 
wi th the effects of Gal actic rotation, and is supported 
by ICarnev et al.l (l2005b who found a significant defi- 
ciency in the number of objects moving in retrograde 
Galactic orbits, at least for metal-poor binary stars in 
the Solar neighborhood. The range of shifts that in- 
crease the correlation amplitude is wide reflecting the 
broad parameter space of velocities and kinematical 
ages represented within the HMXB class. Dynamical 
ejection from the cluster prior to the supernova phase 
can lead to runaway velocities of the order of 150- 
200 km s ' dPoveda et al.lll967l: iGies & Bolton| [l986l 



Pflamm-Altenburg & Kroupal2010l) . This would make 
it difficult to retrace the trajectory of the HMXB back 
to its parent association, and the migration distance 
would be larger than expected from a kick alone. This 
is another factor contributing to the wide range in the 
shifts that increase the amplitude of ^. 

Thus, the observed distribution of HMXBs in the 



Milky Way is consistent with the vie w that these sys- 
tems h ave high velocities, on average. IChevalier & Ilovaiskv 
( 119981) arrived at the same conclusion by deriving the 
pecuHar velocities of 17 HMXBs using HIPP A RCOS 
observations (see also 
concerning unreliable 



van den Heuvel et al 



(l2000l) 

distance estimates for these 
sources). This feature is not unique to the Milky Way 
since HMXBs in the Small Mage llanic Clou d also 
possess large velocities, on average (ICoell2005h . 

Minimum separation distances can provide clues 
to the HMXB-OB co nnection in specific cases (e.g.. 



van den Heuvel et al.l 12000) . For example, there 



are 8 HMXBs in our sample for which the uncer- 
tainty on the line-of-sight distance is smaller than 
the distance separating it from its nearest OB asso- 
ciation: IH 1249-637, Cyg X-3, EXO 2030-H375, 
gam Cas, IGR J18027-2016, IGR J18483-0311, 
IGR J19140+0951, and SS 433. Those with sepa- 
rations smaller than the average minimum separation 
distance (~0.4kpc) are prime examples of HMXBs 
whose migration distances are consistent with those 
expected from post-natal kicks. Small separation dis- 
tances can result from kick velocities that are lower 
than average, which can occur during the formation of 
a black hole since less material is expected to be ex- 
pelled during t he superno va (or gamma-ray burst) that 
created it (e.g.. iMirabel & RodrigueS:2003k) . 

A separation distance larger than 0.4 kpc could im- 
ply that the binary inherited a higher kick velocity than 
average, or that it experienced dynamical ejection prior 
to its supernova (placing the system further away from 
its OB association than would be expected from its 
velocity and kinematical age). Another possibility is 
that the OB association that produced the HMXB is 
no longer active enough to be catalogued: massive 
stars can photo-evaporate their molecular clouds o n 
timescales of 30 Myr (e.g.. lWilliams & McKeil997l) . 

On the other hand, if the HMXB distance is rel- 
atively accurate (i.e., within ±0.5 kpc), and if the 
distance to the nearest OB association is less than 
this uncertainty, then this increases the likelihood 
that the objects were linked in the past. This is the 
case for 12 HMXBs in our sample: 4U 1700-377, 
Cyg X-1, GX 301-2, GX 304-1, H 1145-619, 
1GRJ01583+6713,RXJ0146.9+6121,RXJ0440.9+4431, 
RX J1826.2-1450, SAX J1818.6-1703, VelaX-1, and 
XPer 

Table[3] compiles these results. This assumes (per- 
haps unjustifiably) that the nearest OB association is 
the most likely birthplace of the HMXB being consid- 
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ered. Confirmation of the link (or lack thereof) be- 
tween an HMXB and the nearest OB association is 
only possible through a detailed study of the veloci- 
ties of the objects, which is beyond the scope of this 
paper 

We close this section by pointing out that of the 79 
HMXBs in our sample, 9 are located more than 1 kpc 
from their nearest OB association. Thus, the vast ma- 
jority (89%) of the HMXBs in our sample are situated 
clos e (<1 kpc) to an OB associa tion. This is consistent 
with IZinnecker & Yorkd (l2007b who estimate that 5- 
10% of all OB stars form away from any cluster or as- 
sociation. Distance uncertainties are reported for 6 of 
9 such sources, and these uncertainties exceed 1.5 kpc, 
except for H 0115+634 which has an uncertainty of 
+ 1 kpc. This BEXB is situated more than 2 kpc from 
its nearest OB association. This could represent a case 
where the young stellar bin ary was formed in isolation 
(IZinnecker & Yorkell2007l) . 



6. Summary & Conclusions 

This is the first time that the clustering of HMXBs 
and OB associations in the Milky Way has been 
demonstrated statistically. Since the correlation func- 
tion relates objects in Cartesian space (adding a third 
dimension is trivial), this is a more robust spatial rela- 
tion than those derived in prior studies which focused 
on the distribution of a single dimension such as lon- 
gitudes or galactocentric radii. 

Not only does the correlation function confirm the 
expected view that HMXBs and OB associations are 
clustered together, the characteristic scale of the cor- 
relation contains the vestiges of stellar and Galactic 
evolution. Migration due to the kick velocity gained 
by the binary after the supernova can be constrained to 
a few hundred parsecs. This translates to kinematical 
ages (time spanning the supernova and HMXB phases) 
of ~4 Myr, which is consistent with theoretical expec- 
tations. We point out that these results are based on 
the distributions of the HMXB and OB populations 
which compose the "grand design" of the Milky Way. 
Yet even with this global view, the correlation function 
allows us to deduce the effects of local perturbations 
such as kick velocities. 

The correlation function opens novel areas of re- 
search. Our upcoming Chandra survey of the Norma 
Arm will help locate dozens of new HMXBs which 
will allow us to examine the clustering of HMXBs and 
OB associations within a specific arm. The structure of 



these spiral arms s hould come into sharpe r focus with 
the Gaia mission ( Ferryman etalJ 200lh. The nex t 
hard X-ray surveyor, NuSTAR ( Harrison et al. 2010l) . 
will provide unprecedented sensitivity at these ener- 
gies in the hunt for HMXBs. An increase in the dis- 
covery space of HMXB populations should allow us to 
probe deeper into the evolutionary history of massive 
stars and compact objects. This will permit a better 
understanding of the stellar content and its distribution 
in the Galaxy. 
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Table 1 

Distances reported for the high-mass X-ray binaries used in this analysis 



Mil tnp' 


7 


h 




lA 0535+262 


181.445 


-2.644 


9+0.4 

—0 7 


lA 1118-615 


292.500 


-0.892 


5 ±2 


IE 1145.1-6141 


295.490 


-0.010 


8.5 ± 1.5 


IH 1249-637 


301.958 


-0.203 


0.392 ± 0.055 


3A0114-H650 


125.710 


+2.563 


4.5 ± 1.5 


3A 0656-072 


220.128 


-1.769 


3.9 + 0.1 


3 A 1845-024 


30.421 


-0.406 


10 


3A 2206-H543 


100.603 


-1.106 


2.6 


4U 1036-56 


285.350 


+ 1.431 


5 


4U 1700-377 


347.754 


+2.173 


2.120 + 0.343 


4U 1909+07 


41.896 


-0.811 


7±3 


AX J1820.5-1434 


16.473 


+0.068 


8.2 + 3.5 


Cen X-3 


292.091 


+0.337 


10+ 1 


CygX-1 


71.335 


+3.067 


2.1+0.25 


Cyg X-3 


79.845 


+0.700 


7 9+0.3 

—U.J 


EXO 0331+530 


146.052 


-2.194 


7.5 ± 1.5 


EXO 2030+375 


77.152 


-1.242 


7.1 ±0.2 


gam Cas 


123.577 


-2.148 


0.117 + 0.012 


Ginga 1843+009 


33.038 


+ 1.690 


12.5 + 2.5 


GRO J 1008 -57 


283.000 


-1.822 


5 


GT 0236+610 


135.675 


+ 1.086 


2.5 


GX 301-2 


300.098 


-0.035 


3.5 + 0.5 


GX 304-1 


304.103 


+ 1.247 


2.4 ± 0.5 


HOI 15+634 


125.924 


+ 1.026 


8 ± 1 


H 1145-619 


295.611 


-0.240 


3.1+0.5 


H 1417-624 


313.021 


-1.598 


6 + 5 


H 1538-522 


327.419 


+2.164 


4.5 


H 1907+097 


43.744 


+0.476 


5 


IGR J00370+6122 


121.221 


-1.464 


3 


IGR JO 1363 +66 10 


127.394 


+3.726 


2 


IGR J01583+6713 


129.352 


+5.189 


4 + 0.4 


IGR J06074+2205 


188.385 


+0.814 


1 


IGR J08262-3736 


256.425 


+0.296 


6.1 


IGR J08408-4503 


264.040 


-1.952 


2.7 


IGR Jl 1215-5952 


291.893 


+ 1.073 


6.2 


IGR Jl 1305-6256 


293.945 


-1.485 


3 


IGR Jl 1435-6109 


294.881 


+0.687 


8.6 


IGR J13020-6359 


304.089 


-1.120 


5.5 ± 1.5 


IGR J14331-6112 


314.846 


-0.764 


10 


IGR J 16 195-4945 


333.557 


+0.339 


4.5 ± 2.5 


IGR J 16207 -5 129 


332.459 


-1.050 


6 1+'*'* 



classification 



reference 



BEXB (P, QPO, C, T) 
BEXB (P, QPO, C, T) 
SGXB (P) 
BEXB (P?) 
HMXB (P, E) 
BEXB (P C, T) 
BEXB (P, T) 
SGXB (NS, QPO) 
BEXB (P, T) 
SGXB 
SGXB (P) 
BEXB (P, T) 

SGXB (P, E, QPO, C, VHE) 

SGXB (BHC, muQSO) 

SGXB (BHC, muQSO, QPO) 

BEXB (P, QPO, C, T) 

BEXB (P, QPO, T) 

BEXB 

BEXB (P, T) 

BEXB (P, C, T) 

BEXB (BHC, muQSO, VHE) 

HMXB (P, C, T) 

BEXB (T) 

BEXB (P, QPO, B, C, T) 

BEXB (P, T) 

BEXB (P, T) 

SGXB (P, E, C) 

SGXB (P, C, T) 

SGXB (P) 

BEXB (T) 

BEXB (P, T) 

BEXB 

SGXB 

SFXT 

SFXT (P) 

BEXB (T) 

BEXB (P, T) 

BEXB (P, T) 

HMXB 

SFXT 

SGXB (E?) 



Steele et al. (1998) 
Janot-Pacheco et al. (1981) 
Ray & Chakrabartv (2002) 
Megier et al. (20091 
Farrell et al. (20081 
McBride et al. (2006) 
Grimm et al. (2002) 
Blav et al. (2006) 
Grimm et al. (20021 
Megier et al. (20091 
Morel & Grosdidier (2005) 
Kinugasa et al. (1998) 
Hutchings et al. (1979) 
Ziolkowski (20051 
Ling et al. (20091 
Negueruela etal. (19991 
Wilson et al. (2002) 
Megier et al. (20091 
Israel et al. (20011 
Grimm et al. (20021 
Grimm et al. (20021 
Kaper et al (20061 
Parkes et al. (19801 
Negueruela & Okazaki (2001) 
Stevens et al. (1997) 
Grindlav et al. (19841 
Clark (2004) 
Cox et al. (20051 
Negueruela & Reig (2004) 
Reig et al. (2005b) 
Kaur et al. (2008) 
Masetti et al. (2006c1 
Masetti et al (20101 
Levder et al. (20071 
Masetti et al. (2006b1 
Masetti et al. (2006b} 
Masetti et al. (2009J 
Chernvakova et al. (20051 
Masetti et al (2008) 
SidoUetaL (2005) 
Nespoli et al. (20081 



Table 1 — Continued 



Name 



distance [kpc] 



classification 



reference 



IGR J16318-4848 335.617 -0.448 

IGR J 16320-4751 336.330 +0.169 

IGR J16327-4940 335.093 -1.132 

IGR J16358-4726 337.099 -0.007 

IGR J 16393 -4643 338.002 +0.075 

IGRJ16418-4532 339.189 +0.489 

IGR J 16465 -4507 340.054 +0.135 

IGRJ16479-4514 340.163 -0.124 

IGR J 17252-36 16 351.497 -0.354 

IGR J17391-3021 358.068 +0.445 

IGR J 17404-3655 352.626 -3.273 

IGR J17544-2619 3.236 -0.336 

IGR J 18027 -20 16 9.420 +1.036 

IGR J18214-1318 17.681 +0.485 

IGR J18406-0539 26.670 -0.173 

IGR J18410-0535 26.764 -0.239 

IGR J 18450-0435 28.138 -0.657 

IGRJ18483-0311 29.750 -0.745 

IGR J191 13+1533 49.010 +2.758 

IGR J19140+0951 44.296 -0.469 

IGR J2 1347 +4737 92.172 -3.117 

KS 1947+300 66.099 +2.083 

OAO 1657-415 344.369 +0.319 

PSRB1259-63 304.165 -1.011 

RXJ0146.9+6121 129.541 -0.800 

RXJ0440.9+4431 159.847 -1.270 

RXJ1826.2-1450 16.882 -1.289 

SAX J1818.6-1703 14.080 -0.704 

SAX J2 103. 5 +4545 87.130 -0.685 

SS 433 39.694 -2.245 

SWIFT J2000.6+3210 68.986 +1.134 

VelaX-1 263.058 +3.930 

XPer 163.081 -17.136 

XTEJ1543-568 324.955 -1.475 

XTEJ1810-189 11.362 +0.059 

XTEJ1829-098 21.697 +0.279 

XTEJ1855-026 31.076 -2.092 

XTEJ1946+274 63.207 +1.396 



1.6 
3.5 
2 

7± 1 
10.6 
13 

^-"-5,7 

2.8!g 

8.2.5 

2.7 
9.1 
3.2+1 
12.4 ±0.1 
8 ±2 

I. 1 
3 2+2 

3.6 
2.83 + 0.05 

9.1 
3.6 ± 0.04 

5.8 
9.5 + 1.1 
7.1 ± 1.3 

9 7^+0.61 
^■'-'-0.49 
2.3+0.5 

3.3+0.5 
2.5+0.1 
2.1 +0.1 
6.8+2.3 
5.5+0.2 
8 

1.4 + 0.5 
0.801+0.138 
10 

II. 5 
10 
10 

9.5+2.9 



HMXB (sgB[e]) 
SGXB (?) 
SGXB (T) 
HMXB (P, Be?) 
SGXB (?) 

SGXB (? E, SFXT?) 
SFXT (?, E?) 
SFXT (E, ??) 
SGXB (? E) 
SFXT 
HMXB 
SFXT (NS?) 
SGXB (? E) 
SGXB 

BEXB (muQSO?) 

SFXT (?) 

SFXT 

SFXT (?) 

HMXB (sgB[e], T) 

SGXB (NS?) 

HMXB (T, SFXT?, Be?) 

BEXB (? T) 

SGXB (? E) 

BEXB (? VHE) 

BEXB (? T?) 

BEXB (?) 

SGXB (BHC, muQSO, VHE) 
SFXT 

BEXB (? Q?0, T) 
SGXB (BHC, muQSO) 
BEXB (? T) 
SGXB (? E, C) 
BEXB (? C) 
BEXB (? T) 
HMXB (NS, T) 
HMXB (?, T, Be?) 
SGXB (? E, T) 
BEXB (? C, T) 



Rahoui et al. (2008) 
Rahoui et al . (2008) 
Masetti et al . (2010) 
Liu et al. (2006) 
Chatv et al. (2008) 
Chatv et al. (2008) 
Nespoli et al. (2008) 
Nespoli et al. (2008) 
Mason et al. (2009) 
Rahoui et al. (2008j 
Masetti et al. (2009) 
Rellizza et al. (2006) 
Torreion et al. (2010) 
Butler et al. (2009) 



Masetti et al. r2006a1 
NespoU et al. (20081 
Coe etal. (1996) 
Torreion et al. (2010) 
Masetti et al. (2010) 
Torreion et al. (2010) 
Masetti et al. (2009) 
Tsygankov & Lutovinov (2005) 
Audlev et al. (2006) 
Cordes & Lazio (2002) 
Reig et al. (1997) 
Reig et al. (2005a) 
Casares et al. (2005) 
Torreion et al. (2010) 
Reig et al. (2010) 
Blundell & Bowler (2004) 
Masetti et al. (2008) 
Chevalier & Ilovaiskv (1998) 
Megier et al. (2009) 
Grimm et al. (2002) 
Markwardt et al. (2008) 
Halpern & Gotthelf (2007) 
Grimm et al. (2002) 
Wilson et al. (2003) 



Note. — BEXB=Be X-ray binary; BHC=black hole candidate; C=cyclotron; E=eclipsing; HMXB=high-mass X-ray binary; 
muQSO=microquasar; NS=neutron star; ?=pulsations; Q?0=quasi-periodic oscillations; SFXT=supergiant fast X-ray transient; 
sgB[e]=supergiant B-emission line star; SGXB=supergiant X-ray binary; T=transient; VHE=very high energy 



Table 2 

Kick velocities and kinematical ages reported for a few well-known HMXBs 



name 


velocity [km s ' ] 


kinematical age [Myr] 


migration distance [kpc] 


reference 


4U 1700-377 


75 


2.0 


0.2 


Ankav et al. (2001) 


LS5039 


150 


1.1 


0.2 


,Rib6 et al. (2002) 


LSI+61°303 


27 


1.7 


0.05 


Mirabel et al. (2004) 


VelaX-1 


160 


0.1 


0.02 


Gvaramadze et al. (201 1) 



Table 3 

Minimum separation distances between HMXBs and OB associations in selected cases 



name HMXB distance uncertainty [kpc] distance to nearest OB assoc. [kpc] 



HMXB distance uncertainty < distance to nearest OB assoc. 



IH 1249-637 


0.06 


0.1 


CygX-3 


0.5 


0.7 


EXO 2030+375 


0.2 


0.6 


gam Cas 


0.02 


0.08 


IGR J 18027-20 16 


0.1 


0.8 


IGR J 18483-03 11 


0.05 


0.2 


IGR J19140+0951 


0.04 


0.3 


SS433 


0.2 


0.7 


HMXB distance uncertainty > distance to nearest OB assoc. 


4U 1700-377 


0.5 


0.2 


CygX-1 


0.25 


0.2 


GX 301-2 


0.5 


0.2 


GX 304-1 


0.5 


0.1 


H 1145-619 


0.5 


0.2 


IGR J01583+6713 


0.4 


0.4 


RX J0440.9+4431 


0.5 


0.5 


RX J0146.9+6121 


0.5 


0.4 


RX J1826.2-1450 


0.1 


0.004 


SAX J1818.6-1703 


0.1 


0.1 


VelaX-1 


0.5 


0.03 


XPer 


0.3 


0.1 



